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Abstract 

Background: Deep-sea hydrothermal vents provide patchy, ephemeral habitats for specialized communities of 
animals that depend on chemoautotrophic primary production. Unlike eastern Pacific hydrothermal vents, where 
population structure has been studied at large (thousands of kilometres) and small (hundreds of meters) spatial 
scales, population structure of western Pacific vents has received limited attention. This study addresses the scale at 
which genetic differentiation occurs among populations of a western Pacific vent-restricted gastropod, Ifremeria 
nautilei. 

Results: We used mitochondrial and DNA microsatellite markers to infer patterns of gene flow and population 
subdivision. A nested sampling strategy was employed to compare genetic diversity in discrete patches of Ifremeria 
nautilei separated by a few meters within a single vent field to distances as great as several thousand kilometres 
between back-arc basins that encompass the known range of the species. No genetic subdivisions were detected 
among patches, mounds, or sites within Manus Basin. Although /. nautilei from Lau and North Fiji Basins (-1000 km 
apart) also exhibited no evidence for genetic subdivision, these populations were genetically distinct from the 
Manus Basin population. 

Conclusions: An unknown process that restricts contemporary gene flow isolates the Manus Basin population of 
Ifremeria nautilei from widespread populations that occupy the North Fiji and Lau Basins. A robust understanding 
of the genetic structure of hydrothermal vent populations at multiple spatial scales defines natural conservation 
units and can help minimize loss of genetic diversity in situations where human activities are proposed and 
managed. 



Background 

The spatial scales at which individuals within a popula- 
tion interact and the geographic extent of larval disper- 
sal shape the dynamics of marine populations. Dispersal 
capabilities of some species extend across entire ocean 
basins [1], but larval propagules of many other species 
are retained close to their source [2]. Larval develop- 
ment can impose limits on dispersal. Species that brood 
their offspring (direct development) tend to have more 
restricted distributions than species with long-lived, 
planktonic larvae [3], though exceptions exist [4]. 



* Correspondence: andrew.david.thaler@gmail.com 

1 Marine Laboratory, Nicholas School of the Environment, Duke University, 

135 Duke Marine Lab Rd, Beaufort, NC 28516, USA 

Full list of author information is available at the end of the article 



Species that aggregate in small patches may interact and 
reproduce with other individuals in an area encompass- 
ing a meter or less [5] and larvae that lack broad disper- 
sal potential may recruit to their natal population [6]. A 
sampling scheme that fails to account for the localized 
effects of self-recruiting patches may create an appear- 
ance of panmixia, even if substructure exists among 
patches [7]. 

Species dependent on deep-sea hydrothermal vents are 
restricted to patchy, ephemeral habitats that limit the 
areal extent and occurrence of populations. Hydrother- 
mal vent fields are found on mid-ocean ridges, back-arc 
spreading centres, and submarine volcanoes [8]. Organ- 
isms that thrive at vents are supported by chemoauto- 
trophic microbes that metabolize reduced compounds in 
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the vent effluent [9]. Vent habitats are transient, at tem- 
poral scales ranging from days to hundreds of years 
[10], and constituent species may be subject to frequent 
local extinction and recolonization events [11,12]. Survi- 
val of vent species therefore depends on fast growth, 
rapid reproduction, and dispersal abilities that shape the 
diversity and genetic structure of populations [13,14]. 

At mid-ocean ridges, deep-sea hydrothermal vents are 
distributed along roughly linear axes that may function 
as dispersal corridors [15-17]. Geographic populations 
of hydrothermal vent-dependent species can be panmic- 
tic across the extent of their range (e.g., the shrimp, 
Rimicaris exoculata on the Mid-Atlantic Ridge [18-20]) 
but this is not always the case. 

Evidence for isolation-by-distance in vent species has 
sometimes been ambiguous due to small sample sizes 
and inconsistency in the resolution of various genetic 
markers [21]. Considerable evidence exists for geo- 
graphic subdivision associated with geomorphological 
features that affect different taxa to varying degrees. For 
example, the Easter Microplate is associated with isola- 
tion of northern and southern East Pacific Rise popula- 
tions of mussels, but not of polychaete annelids [22,23]. 
A 2000-m long "habitat gap" across the Equator is 
implicated in the isolation of some East Pacific Rise spe- 
cies and variable impedance of gene flow in other spe- 
cies [23,24]. Similarly, a 350-km long ridge offset, the 
Blanco Transform Fault, isolates Juan de Fuca and 
Gorda ridge limpet populations [25]. The same barrier 
interacts with current regimes and is correlated with 
southward unidirectional gene flow in the vent poly- 
chaete Ridgeia piscesae [17]. Life history and behavioral 
attributes of various taxa result in these differing 
responses to shared dispersal barriers [14]. 

Identification of population structure at various spatial 
scales depends in part on the choice of genetic markers. 
For example, amplified fragment length polymorphisms 
were used to test for fine-scale differentiation among 
discrete patches of the tubeworm, Riftia pachyptila, 
separated by as little as 400 m in a venting area along 
the East Pacific Rise, although sample sizes were small 
(n < 15 per site [26]). More conservative mitochondrial 
and nuclear DNA sequences in R. pachyptila revealed 
panmixia at local scales and isolation-by-distance [27] at 
greater geographical scales [28,29]. 

In contrast to mid-ocean ridge systems, limited atten- 
tion has been afforded to the population structure of 
vent organisms from western Pacific back-arc basins. 
These basins are distributed in a non-linear pattern, 
reflecting the complex tectonic history of the region 
[30]. Hydrothermal vents in western Pacific back-arc 
basins are geographically isolated from vents on the East 
Pacific Rise [31]. Regional isolation of species was 
detected among western pacific vents [30]: the Okinawa 



trough and Izu-Ogasawara Arc have a faunal assemblage 
distinct from that of other western Pacific hydrothermal 
vents, and the faunal composition of the Marianna 
Trough is distinct from that of neighbouring basins [30]. 
Vent species tend to be shared among Manus, North 
Fiji, and Lau Basin, but are distinct from species that 
occur at the Okinawa and Marianna Troughs or the 
Izu-Ogasawara Arc [30]. 

Because back-arc basin hydrothermal systems in the 
western Pacific are located on isolated ridge segments 
(in contrast to the linear, semi-continuous series of seg- 
ments on mid-ocean ridges), it has been hypothesized 
that reduced connectivity among western Pacific back- 
arc basins may yield more endemic vent fauna within 
discrete back-arc basins [32]. Some species endemic to 
these basins appear to be panmictic across multiple 
basins (e.g., the mussel Bathy modiolus brevior [33]), 
whereas others are restricted to single basins (e.g., neo- 
verrucid barnacles [34]). Provannid snail species in the 
genus Alviniconcha represent a cryptic species complex 
composed of at least three evolutionary lineages, one 
that occurs at hydrothermal vents in North Fiji Basin, 
one that is restricted to vents in the Marianna Trough, 
and one that co-occurs in both Manus and North Fiji 
Basin [35]. A similar pattern of strong genetic differen- 
tiation may exist within other species. To date, compre- 
hensive efforts have not been made to characterize 
population structure within vent taxa of western Pacific 
back-arc basins. 

IJremeria nautilei is a provannid gastropod that occurs 
in Manus, North Fiji, and Lau Basins and depends on 
sulphur-oxidizing bacterial endosymbionts for nutrition. 
Sessile adults live in discrete patches near the effluent of 
diffuse-flow hydrothermal vents [36,30]. Females possess 
a specialized brood-pouch in their foot and they release 
ciliated pre-veliger larvae (Waren's larvae) that are 
hypothesized to have long-distance dispersal capabilities 
[37]. Preliminary studies indicated that /. nautilei exhi- 
bits distinct mitochondrial haplotypes in Manus and 
North Fiji Basins [38], but population structure has not 
been assessed at smaller spatial scales-among vent fields 
within basins (henceforth sites), among sulphide mounds 
within vent fields (henceforth mounds), or among dis- 
crete patches on vent mounds (henceforth patches). 

We examined genetic population structure of Ijre- 
meria nautilei from hydrothermal vents in Manus, 
North Fiji, and Lau Basins at multiple scales, ranging 
from meters to thousands of kilometres (Figure 1). A 
nested sampling strategy was employed within Manus 
Basin to test the null hypothesis that /. nautilei exhibits 
no population structure among discrete patches at spa- 
tial scales of meters to 40 kilometres. The entire Manus 
Basin population was then compared to North Fiji and 
Lau basin samples to assess the relationship between 
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Figure 1 Ifremeria nautilei sampling locations in the western 
Pacific. Grey circles are approximate location of sampling. Dashed 
ines represent subduction zones. Three sites sampled in Manus 
Basin are shown in inset. 



increasing spatial scales (1000 kilometres to 3500 kilo- 
metres) and genetic differentiation. Genetic markers for 
differentiation at these scales included partial sequences 
of mitochondrial cytochrome-c-oxidase subunit I, and 
an array of nuclear DNA microsatellite loci [39]. By 
comparing these two types of molecular markers, we 
can separate evolutionary processes, revealed by COI 
sequence data and dependent on mutation rates, from 
ecologic processes, revealed by microsatellite allele fre- 
quencies and based on the recombination of alleles with 
each generation [40]. If the specialized Waren's larvae 
produced by Ifremeria nautilei are adapted for long-dis- 
tance dispersal [37], population structure should be 
minimal over all scales. Alternatively, if /. nautilei dis- 
perse in a manner consistent with other sessile inverte- 
brates with specialized habitat needs [41] genetic 
differentiation may occur at spatial scales less than one 
kilometre. 

Results 

Summary statistics 

Ifremeria nautilei were sampled from Manus, North Fiji, 
and Lau Basins (275 individuals total; Figure 1, Table 1). 
Thirty-six partial COI haplotypes (404 bp) were identi- 
fied from 158 Ifremeria nautilei sampled from Manus 
Basin; an additional 25 haplotypes were identified from 
117 individuals from North Fiji and Lau Basins. Haplo- 
type diversity (H d ) ranged from 0.59 to 0.95 (Table 2) 
and nucleotide diversity (tt) ranged from 0.003 to 0.007. 
Indices of genetic diversity and tests for selection are 
reported in Table 2. Fu's F s values were negative, consis- 
tent with allelic excess driven by recent population 
expansion (but may also be indicative of a selective 
sweep) [42], with two exceptions (Patch 3 and Patch 16). 



Table 1 Ifremeria nautilei sampling locations from Manus, 



North Fiji, and Lau Basin. 


Basin Site 


Mound 


Latitude 


Longitude 


Depth 

(m) 


Manus Solwara 8 


Active 
1 


3° 43.740'S 


151° 40.404'E 


1720 




Active 
2 


3° 43.824'S 


151° 40.458'E 


1710 




Active 
3 


3° 43.668'S 


151° 40.872'E 


1650 


Solwara 1 


Active 
4 


3° 47.436'S 


152° 5.472'E 


1530 




Active 
5 


3° 47.370'S 


152° 5.778'E 


1490 




Active 
6 


3° 47.370'S 


152° 5.616'E 


1480 


South Su 


Active 
7 


3° 48.564'S 


152° 6.144'E 


1300 




Active 
8 


3° 48.492'S 


152° 6.186'E 


1350 




Active 
9 


3° 48.432'S 


152° 6.306'E 


1320 


North White Lady 
Fiji 




16° 

59.950'S 


1 73° 54.950'E 


1971 


White 
Rhino 




16° 

59.950'S 


1 73° 54.950'E 


1971 


Mussel Hill 




16° 

59.950'S 


1 73° 54.950'E 


1971 


Lau Kilo Moana 




20° 

03.230'S 


176° 

08.01 0'W 


2620 


Tu'i Malila 




20° 

59.350'S 


176° 

34.1 00W 


1884 



All sites from North Fiji Basin are within 100 m of each other and marked 
with the same coordinates. 



Nine microsatellite loci were amplified from 40 to 66 
individuals per site within Manus Basin (Additional File 
1, Table SI). Total alleles per locus ranged from 3 to 23 
(mean = 6.6). In permutation tests, allelic richness (Rs) 
did not vary significantly among patches, mounds, or 
sites (10,000 permutations, P > 0.05; Additional File 1, 
Table SI). Three loci {Ifr040, IfrOS2, and Ifr078) were 
monomorphic at the patch level but were polymorphic 
among patches. 

Only eight microsatellite loci were amplified from 20 
to 38 individuals per site from North Fiji and Lau Basin 
(Additional File 1, Table S2). One locus (Ifr086) failed to 
amplify in any North Fiji or Lau Basin samples. The 
total number of alleles per locus ranged from 2 to 23 
(mean = 6.2). In permutation tests, allelic richness (Rs) 
did not vary significantly among sites or basins (10,000 
permutations, P > 0.05, Additional File 1, Table S2). 

Microsatellite marker quality 

Tests for HWE deviation were used to assess the quality 
of sampled microsatellite markers. In Manus Basin, het- 
erozygote deficiency was detected in one locus (Ifr086; 
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Table 2 COI summary statistics for samples of Ifremeria 
nautilei collected from Manus, North Fiji, and Lau Basin, 
and divided by patches, mounds, and sites within Manus 
Basin. 





W 


H 


"d 


F s 


Manus Basin 


158 


36 


0.83 (0.02) 


-27.301 


Solwara 8 


-1-1 


16 


0.80 (0.06) 


-9.719 


Mound 1 


15 


6 


0.74 (0.09) 


-1.208 


Mound 2 


17 


9 


0.83 (0.09) 


-4.337 


Mound 3 


12 


6 


0.82 (0.10) 


-1.272 


Solwara 1 


58 


16 


0.79 (0.04) 


-8.427 


Mound 4 


11 


8 


0.95 (0.05) 


-3.723 


Patch 1 


7 


6 


0.95 (0.10) 


-3.027 


Patch 2 


■-1 


4 


1.00 (0.18) 


n.d. 


Mound 5 


26 


7 


0.71 (0.07) 


-1.377 


Patch 3 


15 


4 


0.66 (0.08) 


0.503 


Patch 4 


7 


6 


0.95 (0.10) 


-0.780 


Patch 5 


-1 


3 


0.83 (0.22) 


n.d. 


Mound 6 


21 


9 


0.80 (0.06) 


-3.403 


Patch 6 


6 


3 


0.60 (0.22) 


n.d. 


Patch 7 


3 


3 


1 .00 (0.27) 


n.d. 


Patch 8 


12 


6 


0.76 (0.12) 


-1.475 


South Su 


56 


20 


0.88 (0.03) 


-4.918 


Mound 7 


18 


11 


0.91 (0.05) 


-5.865 


Patch 9 


6 


5 


0.90 (0.16) 


n.d. 


Patch 10 


2 


2 


1.00 (0.50) 


n.d. 


Patch 11 


2 


2 


1.00 (0.50) 


n.d. 


Patch 12 


8 


6 


0.93 (0.08) 


-2.401 


Mound 8 


22 


12 


0.91 (0.04) 


-5.473 


Patch 13 


5 


5 


1.00 (0.13) 


n.d. 


Patch 14 


10 


7 


0.93 (0.06) 


-2.906 


Patch 15 


7 


3 


0.67 (0.16) 


-2.354 


Mound 9 


16 


8 


0.80 (0.09) 


-2.474 


Patch 16 


7 


-1 


0.81 (0.13) 


1.081 


Patch 17 


9 


7 


0.92 (0.09) 


-2.952 


North Fiji Basin 


81 


20 


0.70 (0.05) 


-14.603 


White Lady 


25 


11 


0.83 (0.06) 


-5.655 


White Rhino 


27 


9 


0.65 (0.10) 


-3.321 


Mussel Hill 


29 


8 


0.61 (0.10) 


-3.033 


Lau Basin 


36 


10 


0.59 (0.09) 


-4.918 



N = number of individuals, H = number of haplotypes, H d = haplotype 
diversity (standard deviation), F s = Fu's F s with significance (P < 0.05) indicated 
in bold. F s are only reported for patches with a sample size greater than or 
equal to 7. 



Additional File 1, Table SI). Heterozygote excess was 
detected in only one locus at the site level (Solwara 1, 
Ifr040; Additional File 1, Table SI). Significant deviation 
from HWE was not detected at any other patch, mound, 
or site from Manus Basin (Additional File 1, Table SI). 

Two microsatellite loci deviated significantly from 
Hardy- Weinberg Equilibrium at sites within North Fiji 
Basin {Ifr068 and Ifrl03; Additional File 1, Table S2). 
Four microsatellite loci were not in equilibrium at the 



basin level (Ifr068, Ifr078, Ifr093, and Ifrl03; Additional 
File 1, Table S2). Neither directional nor balancing 
selection was detected among microsatellites at any spa- 
tial scale within Manus Basin (LOSITAN, P > 0.05), but 
one microsatellite locus (Ifr043) was under positive 
selection (LOSITAN, 25,000 simulations, P < 0.001) at 
the basin level. 

Microsatellite marker identity and excluded markers 

Identity tests were used to assess the utility of each 
microsatellite marker set. Within Manus Basin, probabil- 
ity of identity tests {Pid) an d probability of sibling iden- 
tity tests (Psib) indicated that the nine microsatellite 
markers identify individuals {Pid = 1-5 x 10~ 6 ), including 
those that shared 50% genetic similarity {Psib = 3.5 x 
10~ 3 ). IMa coalescent models require that microsatellite 
markers adhere to the stepwise mutation model; only 
four of nine microsatellite markers {Ifr043, Ifr052, 
Ifr078, and Ifr086) adhered to this model and could be 
used for IMa analysis. Identity tests for these four mar- 
kers suggested that they are insufficient for assessment 
of population structure {P JD = 1.1 x 10' , P S ib = 0.12). 

Within North Fiji and Lau Basins, four of nine micro- 
satellite markers {Ifr043, Ifr068, Ifr086, and Ifrl03) failed 
to amplify, were out of equilibrium, or were under 
selection. These markers were excluded from all ana- 
lyses involving North Fiji and Lau Basins. The five 
remaining microsatellite loci could identify individuals 
{Pid = 4.0 x 10" ), even those that share 50% genetic 
similarity {Psib = 4.4 x 10 2 ). Only two of those micro- 
satellites {Ifr052 and Ifr078) adhered to the stepwise 
mutation model and could be used in IMa analyses. 
Identity tests for these two markers suggested that they 
are insufficient for assessment of population structure 
{Pid = 0.13, P S m = 0.35). 

Population structure within Manus Basin 

Mitochondrial genealogies revealed two frequent haplo- 
types at all Manus Basin sites (individuals per haplotype 
> 30; Figure 2). Less abundant haplotypes radiated from 
the dominant haplotypes in a star-like pattern (many 
shallow branches radiating from numerically dominant 
haplotypes). South Su contained the most private haplo- 
types (k = 10), followed by Solwara 1 (n = 8) and Sol- 
wara 8 (« = 7). 

No indication of population structure was detected at 
the patch, mound, or site level within Manus Basin: 
Analysis of molecular variance (AMOVA) indicated no 
variation among haplotypes using either mitochondrial 
or microsatellite markers {P > 0.05); pairwise compari- 
sons of F S t and 4>st revealed no significant genetic dif- 
ferentiation {P > 0.05); hierarchical analysis of F S t and 
4>st did not reveal any significant population differentia- 
tion among nested samples (HIERFSTAT, P > 0.05; 
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Figure 2 Statistical parsimony network of Ifremeria nautilei haplotypes from Manus, North Fiji, and Lau Basin Area of circles is 
proportional to number of individuals that possess each haplotype. Small black circles represent inferred haplotypes not recovered in this data 
set. Each node represents a one base pair difference between haplotypes. Boxes delineate each putative haplotype group. 



Table 3); assignment tests indicated that all /. nautilei 
collected in Manus Basin constitute a single population 
(Structure, K = 1, Figure 3). 

Population structure among Manus, North Fiji, and Lau 
Basin 

Haplotypes present in Manus, North Fiji, and Lau Basins 
segregated into three groups (Figure 2). Haplotype 
group 1 included all samples from Manus Basin while 
haplotype group 2 contained a mix of individuals from 
North Fiji and Lau Basin (Figure 2). Haplotype group 3 
also contained a mix of individuals from North Fiji and 
Lau Basin but was not directly connected to haplotype 
group 2 (Figure 2). 



Table 3 Hierarchical analysis of F-statistics from 
populations of Ifremeria nautilei sampled within Manus 
Basin using 9 microsatellite loci and a 404-bp region of 
the COI gene sequence and sampled among Manus, 
North Fiji, and Lau Basins using 5 microsatellite loci and 
a 404-bp region of the COI gene sequence. 

Hierachical Level 9 Loci COI 

Patch/Mound 0.03 (0.27) 0.01 (0.17) 

Mound/Site 0.01 (0.77) 0.01 (0.25) 

SiteAotal 0.00 (0.07) 0.01 (0.08) 

Hierachical Level 5 Loci COI 

Site/Basin 0.00 (0.56) 0.01 (0.06) 

Basin/Total 0.05 (0.04) 0.23 (0.03) 

P-values indicated in parentheses. 
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Figure 3 Structure inferred two populations of Ifremeria 
nautilei, finding North Fiji and Lau Basin as one population 
distinct from Manus Basin. Structure inferred two populations of 
Ifremeria nautilei, finding North Fiji and Lau Basin as one population 
distinct from Manus Basin. Manus Basin population indicated in 
black, North Fiji/Lau population indicated in grey. Solid black line 
denotes division between samples from Manus and North Fiji Basin. 
Dashed line indicates division between samples from North Fiji and 
Lau Basin. Eight microsatellite markers were used to test K 
(Additional File 1, Table S2). Likelihood plateaued at K = 2. 



Assignment tests of individual multi-locus genotypes 
identified two geographical regions hosting distinct 
populations of Ifremeria nautilei (STRUCTURE, K = 1; 
Figure 3): one in Manus Basin and a second occupying 
North Fiji and Lau Basins. STRUCTURE output sug- 
gests that one individual from North Fiji Basin may be 
second-generation migrant from Manus Basin and one 
individual from Manus Basin may be a second-genera- 
tion migrant from the North Fiji/Lau Basin population. 
Hierarchical analyses of F-statistics also detected signifi- 
cant differentiation at this regional level (HIERFSTAT, P 
< 0.05; Table 3), but no significant differentiation was 
detected at lower levels [among sites within regions, 
among sites within the Manus Basin, among mounds 
within Manus sites, or among patches within Manus 
mounds (Table 3)]. Pairwise F ST and cp ST values did not 
show a significant increase in differentiation with geo- 
graphical distance among samples separated by less than 
1100 km (Table 4). 

Estimates of migration, effective population size, and 
divergence time 

As reported above, the few microsatellites that adhered 
to the stepwise mutation model were insufficient to ade- 
quately assess gene flow. Only the mitochondrial COI 
data yielded consistent results in the IMa runs. Coales- 
cent estimates of migration rate (IMa) could not rule 
out the possibility of equal and bidirectional migration 



Table 4 Pairwise comparison of Ifremeria nautilei 
populations from Manus, North Fiji, and Lau Basins. 





MB 


NFB 


LB 


Manus Basin 




0.053 


0.055 


North Fiji Basin 


0.514 




0.000 


Lau Basin 


0.530 


-0.005 





Pairwise F ST from 5 microsatellite loci reported above the diagonal, pairwise 
(f> SJ from 404-bp COI sequences reported below the diagonal. Significant 
differentiation after correction for multiple tests (P < 0.05) indicated in bold. 



within Manus Basin populations. Due to the high level 
of connectivity between Manus Basin samples, it is unli- 
kely that estimates of gene flow would converge on the 
most likely solution. The posterior probabilities display 
multiple peaks and gradually increasing probabilities and 
thus results from these analyses should be approached 
with caution. 

Across samples of /. nautilei from all three western 
Pacific basins, coalescent estimates of migration rate 
(IMa) suggest that migration between North Fiji and 
Lau Basin is high (Table 5). Consistent with a single 
panmictic population, m-values could not be con- 
strained between North Fiji and Lau Basin (Additional 
File 2, Figure SI). No evidence for migration between 
Manus and either North Fiji or Lau Basin was detected 
(IMa), suggesting that /. nautilei from Manus Basin are 
isolated from North Fiji and Lau Basin {ml = m2 = 0; 
Table 5). Posterior probability estimates for gene flow 
between North Fiji and Lau Basin were inconsistent, 
due to the fact that Ifremeria nautilei from North Fiji 
and Lau Basin are part of a single, undifferentiated 
population. An interpretation of directional migration 
between these two basins should be approached with 
caution. 

For Manus, North Fiji, and Lau Basin samples of Ifre- 
meria nautilei, estimates of effective population size 
using either microsatellite linkage disequilibrium or coa- 
lescent analysis could not constrain population sizes 
among basins, suggesting that within each basin, effec- 
tive population size is functionally infinite. Estimates of 
splitting time place the oldest divergence between North 
Fiji and Manus Basin, with a relatively recent split 
between North Fiji and Lau Basin (Table 5). 

Discussion 

The complete absence of genetic subdivision in popula- 
tions of Ifremeria nautilei at distances up to 1000 kilo- 
metres suggests this species is able to colonize distant 
vent habitats and that the ciliated Waren's larvae pro- 
duced by /. nautilei are adapted for long-distance dis- 
persal, as hypothesized by Reynolds et al. [37]. Despite 
this dispersal potential, a barrier to gene flow exists 
between Manus and North Fiji/Lau Basin populations 
that are separated by 2500 kilometres. Although mito- 
chondrial COI gene sequences and nuclear microsatellite 
loci are informative at different temporal scales, both 
markers indicated identical patterns of population struc- 
ture in /. nautilei, regardless of spatial scale. 

Local population structure within Manus Basin 

Ifremeria nautilei from the three Manus Basin sites 
belong to a single, panmictic population based on mito- 
chondrial COI gene sequences and nuclear microsatellite 
markers. Although patch sizes were generally small, 
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Table 5 Best fit models for coalescent analysis of Ifremeria nautilei populations among Manus, North Fiji, and Lau 
Basin. 

Population 1 Population 2 0 , B 2 6 a ml m2 r h \x 1 *404 r/fj 2*404 

Manus North Fiji 59.3 23.3 0.1 0.0 0.0 2.09 103,000 345,700 

Manus Lau 66.0 12.6 0.4 0.0 0.0 1.55 79,900 256,400 

Lau North Fiji 66.9 66.9 6.7 6.2 27.2 0.25 12,400 41,400 

0i = effective size of population 2, 0 2 = effective size of population 2, 0 a = effective size of ancestral population, m, = migration into population 1 from 
population 2, m 2 = migration in to population 2 from population 1, x — splitting time between populations. Splitting time is calibrated against two mutations 
rates: Ut = 5 x 10" 8 [21], u 2 = 1.5 x IfJ 8 [76] adjusted for the number of base pairs in the sequence. Posterior probability densities are provided in Additional File 
2, Figure SI. 



several patches (particularly in Solwara 1) had a suffi- 
cient sample size to test and reject the hypothesis that 
self-recruiting patches of /. nautilei might create the 
appearance of panmixia within this sample set. 

Ifremeria nautilei from South Su had the highest 
abundance of private haplotypes and private alleles. 
Under a scenario of colonization with subsequent 
migration, this pattern could suggest that South Su 
might serve as a source population that contributes indi- 
viduals to other sites sampled in Manus Basin. This 
directional gene flow is consistent with the path of the 
St. George's Undercurrent, which enters Manus Basin 
from the southeast and travels northwest, encountering 
South Su first, then Solwara 1 and Solwara 8, before 
merging with the Vitiaz Straight Undercurrent to form 
the New Guinea Coastal Undercurrent [43]. This 
hypothesis could be tested with development of addi- 
tional genetic markers that provide more detailed genea- 
logical information than the microsatellite markers used 
in this study and by sampling and analysis of individuals 
from sites further west that are known to support Ifre- 
meria nautilei. Rapid population expansion, as suggested 
by Fu's Fs and the star-like mitochondrial genealogies, 
could account for the emergence of private haplotypes 
and alleles at each sampled site within Manus Basin. 
Alternatively, the negative Fu's F s values could be a 
result of a selective sweep on the mitochondrial genome. 
Additional sequenced-based nuclear markers would be 
needed to rule out this possibility. 

Basin-scale population structure 

Individuals of Ifremeria nautilei sampled from across 
the known range of the species in the south western 
Pacific could be subdivided genetically into two popula- 
tions, one restricted to Manus Basin and one distributed 
throughout North Fiji and Lau Basins. Under an isola- 
tion-by-distance scenario, genetic differentiation is 
expected to gradually increase with distance [27], but 
distance alone does not appear to create a significant 
barrier to gene flow. North Fiji and Lau Basin samples, 
separated by -1000 kilometres, were undifferentiated. 
Coalescent estimates of gene flow suggest that genetic 
isolation of /. nautilei populations between Manus Basin 



and North Fiji/Lau Basins may have existed, with occa- 
sional migration, for several hundred thousand genera- 
tions (Table 5), but without an understanding of 
mutation rates and average generation times in these 
snails, it is impossible to place these estimates in a geo- 
logical time frame. 

The phylogeographic break between populations from 
Manus and North Fiji/Lau Ifremeria nautilei is striking, 
considering the high degree of mixing within each 
population. Phylogeographic breaks are often associated 
with oceanographic features (e.g., geomorphology, 
hydrology) that form effective dispersal barriers for a 
wide range of taxa (e.g., Cape Hatteras [44], Cape Cod 
[45], Easter Microplate [14]). A related provannid snail 
complex, Alviniconcha spp., occurs throughout south- 
western Pacific basins [30] and is comprised of several 
cryptic species [35]. The observed phylogeographic iso- 
lation of /. nautilei populations thus does not reflect a 
pattern that is shared by other southwestern Pacific vent 
taxa. 

The phylogeographic break between Manus and North 
Fiji/Lau populations of Ifremeria nautilei is not likely 
the result of a colonization event. Colonization would 
result in a founder effect, where the founded population 
contains a subset of alleles from the source population 
[46]. The single mitochondrial haplotype shared between 
the Manus and North Fiji/Lau populations is intermedi- 
ate between dominant haplotypes from the two popula- 
tions and may be a product of incomplete lineage 
sorting between formerly connected populations [47]. 
Haplotype Group 3 (North Fiji/Lau) is more closely 
related to Haplotype Group 1 (Manus Basin; Figure 2) 
and it consists of as many missing haplotypes as it does 
actual haplotypes. This abundance of missing haplo- 
types, relative to other haplotype groups, could be the 
result of inadequate sampling, or it may suggest that 
disproportionately more haplotypes in that lineage have 
gone extinct. Our interpretation is that the two popula- 
tions once existed as a single population spanning 
Manus, North Fiji, and Lau Basins and that this popula- 
tion became isolated through a vicariant process that 
remains to be determined. The presence of potential 
second-generation immigrants in each population 
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suggests that isolation might not be complete between 
the two regions- 
While Ifremeria nautilei occurs throughout Manus, 
North Fiji, and Lau Basins, populations of /. nautilei fol- 
low the trend of greater endemism and limited connec- 
tivity hypothesized for species endemic to back-arc 
basin spreading centres [32]. In this context, it is not 
surprising that a population of /. nautilei is distributed 
through North Fiji and Lau Basin, as these two basins 
share genera and species [30]. Water masses tend to be 
retained within Lau Basin, with some movement of 
northwestward flowing undercurrents from Lau into 
North Fiji Basin (Thurnherr, unpublished data, available 
at http://www.ldeo.columbia.edu/~ant/LAUB-FLEX/). 
This undercurrent movement is consistent with the 
weak signal of directional gene flow from Lau into 
North Fiji Basin (Table 5). Likewise, the barrier to dis- 
persal between North Fiji and Manus Basin that limits 
gene flow between /. nautilei populations is consistent 
with a barrier that was hypothesized to restrict species 
dispersal between these two basins [48]. The barrier 
could be caused by geomorphological obstacles {i.e. the 
Vanuatu Archipelago), the lack of depth overlap 
between sites in Manus Basin and sites in North Fiji 
and Lau Basin, or by a yet to be determined oceano- 
graphic feature. 

Conclusions 

Theoretical and experimental studies suggest that spa- 
tially and temporally unstable environments favor broad 
dispersal capabilities [49,50], which in turn should lead 
to shallow or absent population subdivisions. In 
dynamic systems such as hydrothermal vents, where 
habitat availability is unpredictable, survival depends on 
long-distance dispersal of propagules. No significant 
genetic differentiation was found among samples of wes- 
tern Pacific Ifremeria nautilei at the patch, mound, or 
site levels within the Manus Basin. No differentiation 
was observed between samples of /. nautilei collected 
from North Fiji and Lau Basins, which are separated by 
-1000 kilometres. The Manus Basin population of /. 
nautilei is isolated from that of the North Fiji and Lau 
Basins by an unknown process that limits contemporary 
gene flow. 

Reproductive mode and larval type are often poor pre- 
dictors of population structure in marine environments 
[3,40,41]. Species with broad dispersal potential have 
been reported with high levels of differentiation at spa- 
tial scales of a few kilometres or less [51-53], while spe- 
cies that would otherwise be expected to show fine-scale 
population structure have been reported to show sur- 
prisingly high levels of connectivity throughout their 
geographic range [54,55]. The absence of Ifremeria nau- 
tilei population structure at all but the broadest spatial 



scales is consistent with long-distance dispersal and the 
barrier to gene flow between the North Fiji and Lau 
Basin population and the Manus Basin population is 
likely extrinsic and not related to life history 
characteristics. 

Fine-scale spatial sampling and genetic analysis such 
as that used in this study can inform mitigation and 
best-management practices for mineral extraction at 
deep-sea hydrothermal vents. The Solwara 1 site is tar- 
geted for deep-sea mineral extraction [56]. A robust 
understanding of population genetic structure at multi- 
ple spatial scales can define natural conservation units 
that can be used to minimize loss of genetic diversity 
within and among populations of vent-restricted species 
[57]. For Ifremeria nautilei, high rates of gene flow 
among the sampled Manus Basin sites suggests that the 
Solwara 1 vents are likely to be repopulated from other 
Manus Basin localities, including South Su and Solwara 
8. Monitoring of species recovery and genetic diversity 
as the Solwara 1 population recovers after extraction 
operations cease should add insight into the rates at 
which novel haplotypes and alleles accumulate in this 
species, providing a means to estimate the ages and 
sizes of extant populations. The Manus Basin population 
of /. nautilei comprises a genetically distinct unit that 
should be managed separately from the North Fiji/Lau 
population. 

Methods 

Geographic setting, sample collection, and DNA 
extraction 

Ifremeria nautilei were collected from three hydrother- 
mal vent sites in Manus Basin: Solwara 8, Solwara 1, 
and South Su (Figure 1). One to four patches of /. nau- 
tilei were sampled from three mounds at each site 
(Table 1). Samples were collected during June-July 2008 
with an ST212 trenching ROV modified for biological 
sampling. Foot tissue was preserved in 95% ethanol. 
Additional /. nautilei samples were acquired from a 
cruise that occurred during May-June 2005 with ROV 
Jason II from North Fiji and Lau Basins (Figure 1). Foot 
tissue was stored briefly at -20°C and transferred to 95% 
ethanol prior to DNA extraction. Genomic DNA was 
isolated by Chelex-Proteinase-K extraction as described 
in Thaler et al. [39] and extracted DNA was stored at 4° 
C until amplification. 

COI amplification and analysis 

The mitochondrial COI (404-bp segment) region was 
amplified with the species-specific primers CO/- 3 and 
COI-6 [38] as follows: 10-100 ng DNA template, 10 x 
PCR Buffer (20 mM Tris, pH 8.8; 50 mM KC1; 0.01% 
Triton X-100; 0.02 mg/ml BSA), 2 mM MgCl 2 , 0.2 mM 
dNTP's, 0.5 (iM each primer, and 1 unit Taq polymerase 
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(Bioline: Taunton, MA) in a 20 ul final volume. Reaction 
conditions were as follows: 94°C for 1 minute; 30 cycles 
of 92°C for 40 s, 50°C for 60 s, 72°C for 90 s; final 
extension of 72°C for 5 min. Amplicons were verified on 
1.8% agarose gels. To remove unincorporated nucleo- 
tides, 14 ul of PCR product was incubated with 0.2 ul 
10 x ExoAP buffer (50 mM Bis-Tris, 1 mM MgCl 2 , 0.1 
mM ZnS0 4 ), 0.05 ul Antarctic Phosphatase (New Eng- 
land Biolabs: Ipswich, MA), 0.05 ul Exonuclease I (New 
England Biolabs: Ipswich, MA) at 37°C for 60 min fol- 
lowed by 85°C for 15 min to inactivate the enzymes. Bi- 
directional sequencing reactions were performed using 
the manufacturer's protocol for Big Dye Terminator 
Reaction (Applied Biosystems: Foster City, CA). 
Sequenced PCR product was purified using AMPure 
magnetic bead system (Agencourt: Morrisville, NC) fol- 
lowing manufacturer's protocol, analyzed on an ABI 
3730x1 DNA Analyzer (Applied Biosystems Interna- 
tional), and edited with Sequencher version 4.7 (Gene 
Codes: Ann Arbor, MI). Consensus sequences were 
compared against the NCBI GenBank database to con- 
firm species identity [58] and aligned using the MUS- 
CLE alignment algorithm [59]. A sequence for each 
unique haplotype was deposited in GenBank (North Fiji 
and Lau haplotypes - accession # JQ074110 to 
JQ074134; Manus Basin haplotypes - accession # 
JQ074135 to JQ074170). 

Neighbor-joining phylograms of aligned mitochondrial 
sequences were assembled in MEGA version 4 [60] with 
an Alviniconcha sp. 2 as an outgroup. Statistical-parsi- 
mony networks were assembled in TCS version 1.21 
(default settings; [61]). Arlequin version 3.11, [62] was 
used to estimate haplotype {H), nucleotide diversity (tt), 
Fu's F s , and pairwise </> ST . 

Microsatellite methods 

Nine microsatellite markers (Ifr040, Ifr043, Ifr052, Ifr068, 
Ifr078, Ifr086, Ifr093, Ifr094, and Ifrl03) were amplified 
from Manus, North Fiji, and Lau Basin samples follow- 
ing methods reported in Thaler et al. [39]. To assess 
marker quality, allelic richness and divergence from 
Hardy- Weinberg Equilibrium (HWE) were calculated in 
GENEPOP (version 4.0, [63]). Permutation tests to 
determine significant variation in allelic richness were 
conducted in F-stat (version 2.9.3.2; [64]). Departures 
from HWE toward heterozygote excess or deficiency 
were assessed for each locus using GENEPOP exact 
tests. Loci were screened using LOSITAN to test for the 
potential influence of selection (25,000 simulations; 
[65]). Microsatellite markers that showed deviations 
from HWE expectations or found to be under the influ- 
ence of selection were excluded from subsequent ana- 
lyses. Identity tests (P ID and Psib) were used to indicate 
whether a given set of microsatellites contains sufficient 



information to be useful for assessing population struc- 
ture [66,67]. P ID and P S i B were calculated for all useful 
sets of microsatellite markers (Gimlet version 1.3.3; 
[68]). 

Common statistical methods 

Arlequin version 3.11 was used to conduct analysis of 
molecular variance (AMOVA). HIERFSTAT [69,70] was 
used to assess hierarchical 0 ST and F sx at various nested 
scales from patch to basin. Microsatellite Analyser 
(MSA; [71]) was used to identify significant differentia- 
tion between patches, mounds, sites, and basins. Alpha 
levels were adjusts via Sequential Bonferroni correction 
to account for multiple tests [72]. Structure version 
2.3.3 [73] was used to visualize population structure. 
We used an admixture model with no a priori sample 
data and with sampling locations as prior distributions. 
Analyses were conducted with a 100,000 step burn-in, 
1,000,000 Markov chain Monte Carlo repetitions, and 3 
replicates per level from K = 1 to 12. The most likely K 
was identified by the average maximal value of Ln P(D) 
returned by Structure. The program LDNe [74] was 
used in an attempt to estimate effective population size 
based on linkage-disequilibrium between microsatellite 
loci. 

Isolation with migration 

Migration rate (w), effective population size (9), and 
divergence time (x) between populations were estimated 
using the coalescent-based isolation-with-migration 
model implemented in IMa [75]. All estimates were 
scaled on mutation rates (u) that are unknown for COI 
in /. nautilei, so splitting time was calibrated against 
two hypothetical rates: u x = 5 x 10~ 8 (determined theo- 
retically, see [21]), and u 2 = 1.5 x 10~ 8 (borrowed from 
rates in the gastropod, Littorina littorea, see [76]). IMa 
runs were performed on COI and microsatellite data 
among three sites within Manus Basin and across all 
three basins. Only microsatellites that did not deviate 
from expectations for the stepwise mutation model 
could be used for IMa analyses [77]. 

A series of short (< 2,000,000 steps) IMa runs was 
conducted to optimize model parameters and determine 
the efficient priors for full runs. Prior probabilities and 
heating schemes were established at 0 = 100, m = 100, 
and t = 1.5, in a 40 chain geometric model. Effective 
sample size, autocorrelations, and trend plots were mon- 
itored to evaluate convergence. Three independent runs 
were compared to ensure that marginal posterior distri- 
butions had achieved similar solutions and results were 
averaged across the three runs. Generated trees were 
analysed in L-mode for best fit (default settings) and the 
most likely model was determined using Akaike Infor- 
mation Criterion [78] and 2LRR tests. 
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Additional material 



Additional file 1: Table SI - Summary statistics for nine 
microsatellite loci amplified from populations of Ifremeria nautilei 
within Manus Basin, n = number of individuals, a = number of alleles, 
Rs = allelic richness, P A = number of private alleles, HE = expected 
heterozygosity HO = observed heterozygosity (bold = significant 
deviation from HWE, * = significant heterozygote excess, t = significant 
heterozygote deficiency, - = null amplification or monomorphic 
genotype). Patches that contained only 1 allele not shown. Table S2 - 
Summary statistics for eight microsatellite loci amplified from populations 
of Ifremeria nautilei from Manus, North Fiji, and Lau Basin, n = number of 
individuals, a = number of alleles, Rs = allelic richness, P A = number of 
private alleles, HE = expected heterozygosity HO = observed 
heterozygosity (bold = significant deviation from HWE, * = significant 
heterozygote excess, t = significant heterozygote deficiency, - = null 
amplification or monomorphic genotype). Patches that contained only 1 
allele not shown. 

Additional file 2: Figure SI - Posterior probability densities for 
migration of Ifremeria nautilei between basins in the western 
Pacific, based on mitochondrial COI gene region 
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